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Abstract 


Theory of plasma oscillation is applied to the problem of the scattering of 
~1000 Mc waves in the E layer of the ionosphere. It is assumed that a region 
of the ionosphere having an irregular or turbulent character contains “‘ oscillating 
domains’”’ or ‘‘ coherent units’’ in which free electrons perform synchronous 
oscillations. These domains are shown to be effective in scattering waves of 
frequencies ~(c/v) wp, where wp is the plasma frequency, v the mean thermal 
speed of electrons, and c the velocity of light. A rough estimate indicates that 
the scattered wave may be observable for a proper setup in the experiment, and 


that its angular dependence is (sin}2%)', 9 being the scattering angle. 
1. Introduction 


Several theoretical studies [1] [2] concerning the problem of the galactic radio- 
frequency radiation, as well as of the disturbed solar radiation point towards an 
explanation in terms of plasma oscillations excited by turbulent streams of ionized 
gases. These theories are based on the fact that the plasma oscillation, the organized 
oscillation of free electrons as a whole, interacts with the radiation. It is therefore of 
interest to investigate the possibility that electromagnetic waves are scattered by the 
plasma oscillations. 

An electromagnetic wave is not scattered in passing through a medium in which 
atoms (electrons) are fairly uniformly distributed, even though the medium contains a 
great many atoms which individually produce scattered wavelets. On the other hand, 
fluctuations of density in the medium cause the scattering of light, on account of 
imperfect destructive interferences between the scattered wavelets. An incident wave 
of wave length Z can only be scattered if, in the Fourier expansion of density fluctua- 
tions, waves with spatial periods of the order of Z occur. An attempt to explain the 
scattering of meter waves in the E layer of the ionosphere has been made by Villars 
and Weisskopf [3]. They proposed an idea that there are external causes ‘which con- 
stantly produce large eddies of a certain dimension Z) and that these are divided soon 
into smaller and smaller eddies. The dividing process.ends with that eddy size LZ, at 
which the effect of molecular viscosity is large enough to dissipate the energy of the 
eddy into heat. Consequently, their theory is applicable for incident waves of wave 
lengths intermediate between Z, and L,, that is, meter waves. 
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In the present paper we would like to point out that the scattering of higher 
frequency waves (~10°Mc) is’ also made possible by electrons under the action of 
periodic electric potentials associated with the plasma oscillations. The ionospheric 
turbulence as suggested by Villars and Weisskopf will be certain to cause the plasma 
of the ionosphere to start oscillation [4]. When electric fields are introduced by in- 
complete space charge neutralization due to disturbances, the highly mobile free 
electrons automatically respond to the forces in such a way as to shield out the fields, 
and a systematic oscillation of charged region is set up. The frequency of the orga- 
nized oscillation of the electron gas as a whole is known to be [3] Worg=[Op? +(34T/m)]!/2 
*(2z/4), p being the plasma frequency defined by w*,=4ze? (density of electrons)/m 
and A referring to the wave length. The initial deviation in velocity and density dis- 
tribution of electrons from equilibrium due to the disturbance will build up longitu- 
dinal waves of a given frequency not always restricted to the frequency given by the 
above dispersion relation. In time, however, these will tend to get out of phase with 
each other, and only the organized plasma waves for which the frequency is given by 
the dispersion relation continue to contribute to macroscopic averages such as the 
mean electric potential. It should be noted that in the E layer of the ionosphere w,,,~p 
(for large 4)~10’sec™! is two powers of ten above the mean frequency of collision 
which an electron makes with molecules, so that the plasma oscillates many times 
before the organized motion is lost by close collision of electrons with molecules 
(collisions between the electrons are neglected for the electronic densities encountered 
in the ionosphere). The plasma oscillations will be excited most effectively in the 
turbulent region. Electrons which enter the region in which the plasma oscillations 
are highly excited differ from those outside the region in that they suffer an oscilla- 
tory change in their velocities. We shall see that these electrons scatter light at the 
expence of their extra momenta arising from their oscillatory motions and that the 
scattered wavelets from a set of electrons confined to a common trough of the plasma 


potential (coherent unit) interfere constructively. 
2. Scattering by Plasma Oscillation 


In a plasma an assembly of electrons undergoes an organized oscillation, the plasma 
oscillation, brought about by the Coulomb forces between charges. The plasma oscil- 
lation or the plasma wave, can be analyzed into a surperposition of plane waves which 
can have all wave lengths, down to limiting waves whose wave lengths are of the 
same order of magnitude of the Debye length [4] [6]. 

Each individual electron is in a state of continual fluctuation in its motion under 
the influence of a periodic potential associated with the plasma wave, and the elect- 
ron can be deprived of its extra momentum and energy by an electromagnetic wave 
which traverses the medium. That is, a part of the monentum and energy of the 
plasma wave can be transferred into incoming electromagnetic wave through the 


intermediary of electrons. We might describe this by saying that the light is scatter- 


ed by a plasma oscillation. As the plasma is a system of a set of individual electrons 
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interacting with the collective plasma oscillations with the uniform background of 
positive charge, let us imagine as the smallest copy of a plasma an ‘atom’ consisting 
of an electron plus a plasma wave. The plasma is now looked upon as containing a 
great many such atoms which individually can produce scattered wavelets. The energy 
of the atom is degenerate in the sense that a given value of the total energy is dis- 
tributed in different ways on the electron and the plasma wave, although these states 
will be different in momentum from each other corresponding to different distributions 
of the total energy. The scattering process consists in the absorption of a primary 
radiation and the emission of a secondary radiation. The scattering atom may be left 
either in its initial state, or in some other state. In particular, if the atom returns to 
a state which is equal to the initial state in energy but different in momentum, the 
emerging radiation carries away a definite amount of momentum (but not energy) 
from the atom, so that the process will appear as a reflection from the atom. 


3. Determination of the Optimum Frequency 


Since we are interested in the scattering of cm waves whose wave lengths are 
much shorter than the linear dimension of the turbulence such as the eddy, we can 
get a better idea of how the plasma oscillation works in the scattering by taking an 
idealized plasma wave of infinitely long wave length. The plasma wave now varies 
only in time and therefore exchanges energy, but not momentum with electrons. An 
electron acted on by a periodic perturbation of the plasma potential exchanges, say 
absorbs, quantum fw, and makes a transition to a virtual state, in which energy and 
momentum are not related by a free particle equation E=p*?/2m, and backs again to 
the initial state, reemitting Aw,, in a time of the order of t=1/wy. If in the time r 
available from the uncertainty principle a radiation fails upon the electron in the 
virtual state and gives a kick to it so that the electron receives such an amount of 
momentum Adk as to satisfy E+hwyp=(p+hdk)/2m, the free particle equation, the 
electron becomes a real one and stays in the new energy state indefinitely while one - 
quantum is lost from the plasma oscillation, the ‘atom’ making an internal conver- 
sion. The radiation then gets momentum —A4k, and is deflected with an angle of 
deflection, say 9. @ is determined from the above free particle relation, or simply from 


dkv~ +, (1) 
to be 


Ak=2ksin 6 = w»/vcos(dkv), (2) 


where v=p/m, + corresponds to the absorption (+) and the emission (—) of the 
longitudinal quantum, and k refers to the wave number of the incident wave. 

If space occupied by electrons were unlimited, the momentum change 4k of 
electrons would be completely determined by the requirement that it be equal to the 
momentum change of the incident light. However, since we are ‘considering those 
electrons that are contained in a certain region of space of, say, linear size Z (which will 
be equated to the wave length of the plasma wave), 4k is not defined to within 
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(2z/l). Now it may be seen that the solutions 4k, of (1) center densely around the 
value 4k=w,/v with Mee, Therefore setting cos(dky)=1 in (2) gives the angular 


frequency at which the deflection occurs as 


Wp c 


= (3) 
2 sin 
2 
For the E layer of the ionosphere where ,~10’sec! corresponding to the mean 
electron density and v~10’cm sec™!, the mean thermal speed, we have 


w/22z~10°Mc/sec, (4) 


which is higher than the penetration frequency, w,/2z, by a factor (¢/v). 

Lastly we notice from (1) with fixed 4k that if, for instance, electrons of velo- 
city v scatter light by damping the plasma oscillation, those of velocity —v do the same 
by exciting the oscillation. 


4. Interaction Hamiltonian 


The scattering of a light quantum is the second order action of the vector 
potential A of the electromagnetic field. Since w),<w, we start with the interaction 
Hamiltonian [5] 

e 
ciet tS » A*(x:) (5) 
é 2£mc 
which describes the scattering of a light quantum on a set of free electrons in the 
nonrelativistic limiting case. So long as we do not know from which electron the 
light quantum is scattered, we sum over the electron coordinates. ‘ The Fourier compo- 
nent of A(x;) responsible for the scattering k,—k is 


Agi tk koxe _ 4 rhc? i(ko —k)xi (6) 


———(@€) € , 
where @, e€ represent unit vectors in the direction of polarization of two light quanta 
of frequency w and of wave vectors ko, k. 

Our system is a set of individual electrons interacting with the collective plasma 
oscillations. Since a photon, striking an electron and being scattered, causes the 
internal conversion of the atom, that is, a non-reversing absorption or a real emission 
by the electron of the quantum Aw, which must be made up by a change of the 
energy or the momentum of the electron, we consider the expression 


e Vny 


: sre eae [de — i Aki sxi 1 (7) 
Zz 8 


where |k—ky|=p/<v>,v, Vn, is the number of atoms contained in a volume V=2 
that is, the number of individual electron degrees of freedom (+3) and x denotes ne 
summation over all the solutions 4k;, of (1) with v=v;. In case of large Vn, (7) 


reduces to 


Vinge Va x 
i= e Aen dks (8) 
é 
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the summation over 7 being this time only for those Vz,. electrons for which 
| ky) —k— 4k; |< (27/0), (9) 


where 
Ak; ~Op/V;- (10) 

When /->0o (9) simply expresses the law of conservation of momentum. In obtaining 
(8) we have neglected a number of terms which are multiplied by a phase factor with 
a nonvanishing argument (>(2z/l)x), exp[é(k)—k—4k)x,]. There is a strong tend- 
ency for various such terms entering into (7) to cancel out, because atoms (electrons) 
are distributed over a wide variety of positions in the region of linear size 7? and there 
would be just as many atoms in the region where exp[i(k)—k—4k)x;] is positive 
as where it is negative [6]. These argument, however, will not apply when the ex- 
ponential function does not vary appreciably inside the region; hence, the condition 
(9). In this way H. becomes practically independent of the electron coordinates, and 
despite the random location of the electrons the scattered wavelets from various 
electrons all add up cumulatively (interfere constructively) to produce a scattered 
wave. 

One can obtain a simple physical picture of this. An incident wave ¢ which 
falls upon each atom varies with the position of the atom as exp(zkyx;). The scat- 
tering of the light by our atom is like the first order action of a perturbing potential 
U varying in space as exp(—7z4kx,): 


U(xi) = (€?/m)(2 zh /wp)(€p€) exp (—74kx:;) 


The scattered wave is the sum of a set of wavelets originating at the positions of the 
atoms. Since each atom contributes according to the product of the potential U(x;,), 
and the incident wave ¢(x,), an outgoing wavelet that spreads out from x; as 
exp[zk(x—x,)] is weighted with the amplitude factor U(x,)¢(x,). Thus the scattered 
wave at a point x is obtained by adding up the contributions of each atom with a 
phase of exp(zkx) exp [2(k)—k— 4k)x;,]. 

Hence, if ky)—k—4k~0 is fulfilled various wavelets from various atoms add in 
phase with each other: 


nA exp (¢kx ) exp[? (Ky) —k— 4k) x,]= Vn, A exp (¢kx) (12) 


We observe that the square of H and hence the scattering cross section is pro- 
portional to m®,,.. Hence our process will bring about a noticeable effect in case of 
large m,.. It should be noted, however, that, since the collective components of the 
density fluctuations of electrons have been expressed in terms of plasma waves, the 
number per unit volume of individual electron degrees of freedom is reduced from 3n 
to 3m,=3n—nyp. Here n is the density of electrons and v, the number of collective 
degrees of freedom given by [6] 


Mp=(41/3)(1/Ap), 


where 4p denotes ‘Debye length.’ Since 4»~1cm in our case, we are justified in 
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setting m,;~n~10'cm-*, We also note that not all individual electrons are effective 
in deflecting incident waves in a given direction, instead a fraction of them, 2,,=2,;42 
/4n, where 42 is a small soltd angle centering about 4dk~(k,)—k) is effective. It 
may be seen that 42/4x~1/AkI,* where 1=V"/3 is the linear dimension of the effec- 


tive volume V in which coherent scattering can take place. 
5. Efficiency for the Scattering 


From H, we obtain in the usual way, the scattering cross section per unit 
volume for the scattering of light by a set of ‘atoms’ 


o/V~rin®,.V=7)2(n/4k)L, (14) 


where a factor of order unity is omitted. 7) is the classical electron radius defined 
by 7»>=e?/mce®. Since the plasma oscillations cannot be sustained for wave lengths 
shorter than a critical distance known as the Debye length, 4p~<v>,4v/wp (which is of 
the order of the distance traveled during the period of an oscillation by an electron 
moving with the mean thermal speed), the actual wave lengths will be much larger 
than Ap. We infer roughly: /=Z,~10°A,, that is, the dimension of the smallest 
eddy. Inserting this into (14) and setting 4k=w,/v~1lcm™! corresponding to the scat- 
tering of 1000 Mc wave give 
o/V~10-"cm"!. (15) 


The useful scattering volume depends upon the geometry for the scattering which 
corresponds to the set up in the experiment. By using the same geometry as Bailey 
et al. (for the scattering of meter waves) the ratio P,/P, of the ‘power received to the 
power transmitted becomes 


P, 40 __bA :  bA 
P, V sin(6/2)D? Shia “ oe D*sin (6/2)’ gf 


where 6b is the thickness of the scattering layer having an irregular or turbulent 
character, A the aperture of the receiving antenna and D the ray length which is 
practically identical with the surface transmission distance. The following set of 
values of the parameters are taken from the paper by Villars and Weisskopf [3]: 
A=3x10%m*, D=10'cm, b=5%x10%cm, 
6=scattering angle=24°. 

Substituting these values, we get 

Pr 19-18 (17) 


t - 
It is of interest to compare (17) with the experimental data for the scattering of 50Mc 
wave which indicate a ratio [3] [7]: 0.36 x 107". ; 
We actually should expect, however, a somewhat = s result than (17), be 
cause the scattering volume will not be filled everywhere with the oscillating domains. 


* From (1) (Ak+2x/l)vcosa=wp=Akv(a=4Q'*) and cosa=Ak/(Ak+2x/l)=1—25/Aki. Since 
cosa=1—«?/2, we get AQ =4n/Akl. 


~) 
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6. Summary 


Our argument was based on the fact that the energy of the ‘atom’ is degene- 
rate with regard to changes of the momentum when the total energy of the system 
of electron plus plasma wave is kept fixed. Under the strong periodic perturbation, 
that is, when the amplitude of the plasma wave becomes large in the atom, reversing 
changes in the distribution of the energy (exchanges of hw, between the electron and 
the plasma wave) occur, so that the electron in the atom must be thought of as 
covering several virtual states centering about a real state at once. We have thus 
far assumed that when a light quantum falles upon the atom a non-reversing 
transition between degenerate levels, which we called the internal conversion, occurs ; 
the momentum thereby liberated is just as large as to deflect the 1000 Mc waves. A 
further possibility to be taken into consideration is that the atom may be left after 
the scattcring in its initial state, when no permanent change in the distribution of 
energy between the electron and the plasma wave (net exchange of hw ,) does occur. 
In other words, the electron in the atom may alone scatter light, the plasma oscilla- 
tion playing no role. Then we must put k,—k~4k=0 in (8) in case of uniform 
distribution of electrons, which means that the light after the scattering travels off in 
the same direction as the incident one. The only effect of the interaction Hamiltonian 
H is therefore that of refraction, a change in the phase velocity as it passes through 
the medium. On the other hand, when there are such fluctuations of the electron 
concentration that there are more electrons in the region where exp[z(k)—k)x,] is 
positive (negative) than where it is negative (positive), H produces a scattered wave 
due to imperfect destructive interferences between the scattered wavelets. Thus the 
presence of density fluctuations whose spatial periods are of the order LZ is essential 
for the scattering of waves of wave lengths ~L<c/w,. This is the content of the 
scattering of meter waves. 

One now obtains a general feature of the scattering in the FE layer of the 
ionosphere of waves above the penetration frequency as follows. An electromagnetic 
wave which enters a region of the ionosphere, in which the electron concentration 
fluctuates, experiences a changing index of refraction and is scattered. When the 
wave length of the incident wave becomes comparable or smaller than L,~10°cm, the 
linear dimension of the smallest eddy, the scattering will fall off rapidly. As the wave 
length decreases further from this value, however, and the frequency approaches three 
powers of ten above the penetration frequency, the wave will again begin to be 
scattered by the plasma oscillations. Thus, tied up with the theory of the scattering 
of meter waves, our argument leads to the selective reflection of electromagnetic 
waves impinging on the ionosphere. 
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Abstract 


The differential equations of the dynamo theory are solved by numerical 
integration, taking into account of the daily variation of the anisotropic conduc- 
tivity. It is shown that the observed wind in the E region is sufficient for the 
requirement of the dynamo theory. The drift of ionization is discussed and it is 
shown that at night some part of electrons in the lower F region will descend 


into the E region in the middle latitudes. 
1. Introduction 


It is widely believed that the small periodic variations in the earth’s magentic 
field are caused by electric current set up by the dynamo action resulting from the 
solar, S, and lunar, Z, tidal motions in the ionosphere ({1], [2]). 

Cowling [3] and Bates and Massey [4] examined the electrical conductivity of the 
ionosphere and it was shown that if the direct electrical conductivity transverse to the 
magnetic field is accepted to be effective to the S or Z current, then the conductivity 
1s much smaller than that required by the dynamo theory. 

Hirono [5] examined the influence of the Hall current to the electrical conduc- 
tivity and the drift of ionization in the ionosphere and showed a possible explanation 
of the electro-jet on the magnetic equator. He showed also that the effective conduc- 
tivity of the ionosphere will be sufficient for the S and Z current on the basis of the 
anisotropic conductivity. The same problem was studied independently by K. Maeda 
[6], Baker and Martyn [7], Cowling [8], Lucas [9], and Fejer [10], and their results are 
essentially similar. 

Recently K. Maeda [11] introduced the daily variation of the anisotropic conduc- 
tivity into the dynamo theory and discussed the drift motions of ionization under an 
assumption of the wind system derived from particular spherical harmonics. Accord- 
ing to the results of those investigations, it is very promising that a number of 
anomalies of the ionosphere is due to the drift motion closely related with the S cur- 
rent system and an exhaustive investigation of the problem is very necessary, although 
at present there are many uncertainties concerning basic physical quantity in the 
region. 

In the present paper the problem is attacked by a fairly general mathematical 
method relating to the previous work on the equator [25] and discussed in some detail, 
based on an atmospheric model mentioned in section 3. 
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2. Differential Equations of the Dynamo Theory 


Consider a right handed cartesian coordinate x, y, 2 with the x axis pointing 
south, the y axis east and the z-axis vertically upward. When the electric field (E,, E,) 
is independent of height in the ionosphere, the horizontal electric current (J,, I,) inte- 
grated with height is given by ([5b], [12]), 


[,=4,E,+4,yEy } (1) 
i, = dy,E,+ayEy 
where 
fic (o,.dh, O,y= (o,,dh, ly = \owdh, Ay. = \vedh. 


. ° 


From (1) we get the electric field (E,, Ey) as follows 


1 1 
SS eS 
E. K, m Ky 
| 1 e 
E,==—L+—=I, 
y Ky a y 


where 
1/K,=ay/4, 1/Kyy= —1/Kyz= @zy/4, } ( 2’) 


1/K,=a,/4, 4=a,ay+a?,,y. 


The values of Z, J, and K, etc. can be evaluated using observed data as shown in 
sections 3 and 4, therefore we can determine (E,, Ey) from (2). The electric field (E,, 
E,) consists of two parts: an electrostatic field (X, Y) and a dynamo field (vH,, —uH,), 


thus 
E,=X+0H, | (3) 


E, = Y— uH, 
where (uz, v) denotes the horizontal air velocity and H, the vertical component of the 
earth’s magnetic field. The electrostatic field (X, Y) is derived from a potential S, 


OS OSs 
en Sa, 4 
ado’ § asin 602 s 


where a stands for the radius of the earth up to the current sheet, @ for the colatitude 
and A for the eastward longitude measured from midnight and, at the same time, the 
local time in angle. It is shown in [5b] that (X, Y) is approximately constant between 
the E and F regions. The horizontal air velocity v=(u, v) can be represented [13] by 


v=u(S)+0(T) (5) 
where curl o(S)=0, div o(T)=0. 


We shall call v(S) and v(T), the scaloidal and the toroidal velocities respectively. It 
would be best to treat our problem using the air velocity represented by (5), however, 
we shall use v(S) or o(T) separately for mathematical simplicity in the following. It 
is assumed that these velocities are constant between 90 and 160km and satisfy the 
condition (33) in higher region. 


(a) Scaloidal Wind System 
In this case the wind velocity is irrotational and can be derived from a scalar 
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potential as in the usual dynamo: theory [14]. 


We take ae Ov sare OP . (6) 
aoog asindoa 


Substituting (4) and (6) into (3) and after some differentiation with @ and 4, we can 
eliminate ¢ from (3) and obtain the following differential equation with S, 


a ee ee ee eee a = 
062 * sin@cos0 00 * sin?@ 02? 00 ' sin@cos@ ' sin@ @A J’ 


The right hand side of this equation can be calculated using the value of (E,, Ey) 
determined by (2), and expressed by 


d3{p,.(8) sin 24 +q,(8) cos nd}. (8) 
If we put S=)>}{f,(9) sin 22 + g,(8@) cosnd}, (9) 


then (7) can be transformed to the following simultaneous equations 


Bp ae 
= in ee Sey, 
dé? sin20 dé sinto”” p,(?) 
d*g,, 2. fia 
dé? sin26 dé sin20 =" 


(10) 


=q,(9) 


where ~, and gq, are known functions. If we solve these equations separately, the 
value of S is obtained up to the necessary terms and from (3) and (4) we get the 
electrostatic field and the wind velocity (zu, v). 
(b) Toroidal Wind System 

In this case the wind velocity is solenoidal and according to Elsasser [15], we 
have vo(T)=(1/a) curl (rg), where r denotes position vector drawn from the centre of 
the earth. Then we have 


a 
“= asind OA 
(11) 
peepee i 
aoe 


Substituting (4) and (11) into (3), and after some differentiation by 6 and-A we can 
eliminate S from (3) and obtain the following differential equation 


0g afdE, 1 OE 2 
Fe cae 0-Ey— = ae sass 
aA Al Ht Bedouin srr rs = } (c =} (12) 


The right hand side of this equation can be determined from (2) and expressed by 

p» {p,(9) sin 24 + q,,(0) cos na}. (13) 
If we put 9= Xi f,(0)sin nd + g,(8)cosna}, (14) 
then (12) is transformed to the following simultaneous equations 


ae 


From these equations, we can determine the value of ¢ up to the necessary terms, and 


(15) 
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from (11) and (3) the distribution of wind and electrostatic field. 


3. Model of the Upper Atmosphere 


At present the structure of the region above 100km level is still uncertain. 
Ferraro [16] discussed the diffasion in the F2 region and concluded that the gas den- 
sity at 300km level is not much léss than 10!°/c.c. Spitzer [17] estimated the tempera- 
ture of the upper ionosphere to be 1500°K from the escape of the helium in the 
region. We take a model M (1500°K; 200km) out of models by Bates [18], in which 
the temperature above 200km is 1500°K, and the gas density at 300km is 10!°/c.c. 
For this model the semi-diurnal variation of the F2 region can be interpreted by 
vertical drift theory [12]. 

Recently Ratcliffe [19] discussed the diffusion under gravity in the F2 region 
and has given strong reason that the gas density at 400km level is not much less 
than 10°/c.c. According to Inoue [20], a steep temperature rise above the E region 
may be due to strong absorption of the solar radiation at the time of dissociation of 
molecular oxygen which has a fairly great concentration above the E region [21]. 
Therefore this model of the atmosphere seems to be better than the low density model 
determined by Havens et al [22]. For the E region which has the maximum electron 
density 1.5x10°/c.c., the density of negative ions being taken as effectively zero [4], 


we have a= { .aidh = 6.85 x 10-9 


(16) 
a,= { adh =2.19x10-8 (e.m.u.) 


where the integrations are carried out between 90 and 160km. For the F region which 
has the maximum electron density 10®/c.c. we have 
B= { adh =2x 10-8 (17) 


e 


where the integration is carried out between 160 and 300km, and \ oodh is much 
smaller than this value. At first sight the electrical conductivity of the F region is of 
comparable order of magnitude to that of the E region. But induction drag will be 
effective in the F region as pointed out by Cowling [3] and later by Hirono [5c] and 
Martyn [7]. Therefore the effective conductivity of the F region will be much smaller 
than that of the E region for the S or Z current system, however, this relation may 
not hold for the current system with shorter period or other modes. At present we 
neglect the conductivity of the F region, i.e. we assumed (33) in section 7, to obtain 
the first approximate solution. The variation with latitude of o,, etc. is investigated 
in detail [5] [6] [7] [9] and Fejer [10] gave a convenient formula as a first approxima- 
tion of (o,.dh etc. with latitude. When the electron density is independent of latitude 
on any level, we denote a, etc. with suffix 1, then according to Fejer [10] it follows 
that 

Gz, = 90,/L, 

Oy; = (aa Sin’*¢ + a,03 cos’¢)/L, (18) 

Azy, = 4,0) Sin $/L, 
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where L=a,cos’¢+a,sin’¢¢, and ¢ denotes geomagnetic dip, 


9 


ia \oudh and a,= \oodh = (0, + 7, )ab 


Numerical calculation gives 
@3=1.52x1077 (e.m.u,), (19) 


where the integration is carried out between 90 and 160km. The values (16) and (19) 
refer to the electron density on the sub-solar point at sun-spot minimum period, and 
in the following we inquire the variation with solar zenith angle. We put the electron 
density N=vN, at any level and N=WN, on the sub-solar point i.e.0<v<1. Then we 
have 


Oz =VOgi5 Ay=VAy, Oey =VEzy1 = — Aye (20) 


Substituting (18) and (20) into (2’) and considering a)>a,, a, we have the following 
values as a fairly good approximation 


-t 1 sin’¢ 
x .. v Ki : 
1 1 sing 
=— . 21 
igs pele (21) 
ete ee eo & 
Ky # Ky 


where K,=a,+a,2/a, and Ky=a,+a,?/a,. Inserting the values of (16), we have 
K,=7.69 10-8, Ky=2.4x 1078. (22) 


According to a more precise calculation the first and second relations of (21) hold on 


Lt : ; 4 1 

the equator, but there we have : , K 
when one goes several degrees from the equator. As (19) and (22) give at iw 1 
a3 


may have somewhat larger value of EF, from (2) and (21). But since we iaas 
the current system (J,, f,) by a few spherical harmonics (see section 4), J, is about 
one half of the observed value in the vicinity of the magnetic equator, and hence we 
shall be able to obtain a fairly good arr field from (2) and (21). Finally we shall 
give the expression of v. We take y= + (ay+a,cos x +a, cos*z), where x denotes zenith 
distance of the sun at equinox and «c= wed, +a,. After a similar calculation to that 
of Chapman [2], we get 
— giao +334. ()cosni. (23) 

The ratio of the electron density at noon to that at midnight is taken to be about 12 
on the equator, then we get the following values: ; 


4=0.6, @,=1.2, a,=0.8 and e=2.6; 
Ay=5.74—0.56sin#, A,=—0.54—0.12sin60 
A, =2.38+0.42sin@, A,=—0.53+0.12sin0 (24) 
A,=0.40—0.10sin 8, 
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4. S, Curren.System at Sun-Spot Minimum Period 
According to Hasegawa and Ota [23], the S, current (J,, I,) at sun-spot mini- 
mum period is derived from the current function YR", where 
Ry =a{snsin mi+encos ma} P™ (0) (25) 


and A denotes the local time in ceverhe: and the values of s™, c™ are shown in table 1 
and P%, denotes associated Legendre function of Schmidt. 


Table 1 
= 2<105° c,!= —9.3x10-* 
$,2= —1.3X10-6 c,2=5xX 10-8 
s3=0.9x10-6 c,8=—-1.3x10-8 


These values were determined from an average of the data of northern summer and 
winter in the Second Polar Year (1932-33). The data analysed came from 46 stations 
located between latitudes +60°, therefore (25) gives a fairly good approximation in the 
range of latitudes. In the polar regions there is yet no determination of mathematical 
expression of the S, field. Moreover the complete current function should be the sum 
of (25) and zonal terms which are not yet determined [24]. For our calculation both 
these unknown expressions should be determined by any means, and we did in the 
following way. Suppose that (J,, J,) is derived from a current function R given by 


R=R(0)+0 3 Rete (26) 


where Rwm,, is defined by (25) with table 1, and R,(@), 8 are functions of 0 only. Then 
we have 


. ; 
has = Pee =)> (Siva 14COS mi —cm, sin mi) EPR) 
asin00A #21 = a 
I,= eae 3 (sR.,sin mA + Ge,,COS md) ——_ 0(0 Pm+1) 
aoe a0 
where 2,(8)= — oR. For mathematical convenience, and under the supposition that 


the current density due to the solar wave radiation would be minimum at the poles, 
we take the following values of 0: 

(1) Polar regions: 9=sin’36 for 0<@<25° and 155°<@<180° 

(II) Moderate and low latitudes: @=1 for 35°<8@<145° 

(III) Boundary region: For 25°<@<35° and 145°<8@<155°, 
the values of @ are determined so as to connect the curves 6 =6(@) in (I) and (II) 
regions continuously up to the derivatives of the second degree with 0, using Bessel’s 
interpolation formula. me 

Suppose that the wind velocity (w, v) is periodic and does not contain any com- 
ponent independent of local time. From the theory of potential, the electrostatic field 
Y should be also periodic for local time, but X is not constrained in this way. Obser- 
vation of (3), then, shows that E, should be periodic but not with E,. Substituting 
(21), (23) and (27) into (2), we can obtain the expressions of E, and £, in Fourier 
series for 4. By equating the zonal term (the term independent of local time) of E, to 
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zero, we can determine the value 4(0) of (27) with colatitude. Introducing this (4) 
into the expressions of EF, and E,, we obtain 
3 3 
E,=Eo,+>){a,(9) sin 2A +8,(9) cos na} = Bort Enc, A) 
*%=1 a= 


3 8 (28) 
Ey= d}{c,(0) sin nd +d,() cos 2a} = EnV, 4) 
n=1 n= 


where a,(9), E,,,(0, 4) etc. are obtained in analytical forms as functions of 6 and/or A, but 
the whole expressions are omitted for brevity. The amplitudes of E,,, E,, and the 
local time when the components attain the positive maximum value, are shown in 
Fig. 1. The value of E,, is shown in 
Fig. 2, and this field is shown to be 
electrostatic from the above arguments. 


3 
¢€ + 50x10" 
= 


a 10 20 30 40 50 60 70 80 90 
CO-LATITUDE 


Fig. 2 The variation of E,, with colatitude 
(independent of local time). 


0 10 20 30 40 §0 60 70 eo = 80 
CO-LATITUDE 


Fig. 1 The amplitudes and phases of E,, and 
Eny. The phases denote the local time when the 
maximum values are attained. 


5. Numerical Method of Solution 


Solutions of the dynamo equations are obtained separately for the two cases in 
section 2. 

(a) Scalotdal Wind System 

In section 2, equations (10) are obtained for the potential functions f, and om 
Substitute (28) into the right hand side of (7) and the values of p, and q, in (8) are 
calculated, for which differentiations with 6 are all carried out by analytical mcthods. 
Thus the right hand side of (10) are obtained. The boundary conditions satisfied b 
f, and g, are: f, and g, must both vanish at the poles (for 6=0 and @=n), fe 


and Je must both vanish at the equator (for 0=2/2). An infinite number of solu- 


tions was obtained, however, for these boundary conditions (see section 6), since (10) 
involves singular points at @=0 and 2/2. But since at the Yc saad H,=0, therefore, 
from (3) we have Y=E, and it follows that 


ie 
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Sn=—(a/n)d,, Ln=(a/n)c,, (29) 
where c, and d, are given by (28). If we add the condition (29) (at 0=7n/2) to the 


above boundary conditions, the solutions of (10) are uniquely determined in the follow- 
ing way. The well known Runge-Kutta’s method is adopted and the integration is 
performed point by point, the range being taken 1°<0<89° for practical purposes. 
First we get a particular integral y, of one of equations (10) such that y,=0 and 
dy,/d0@=k, at @=1°, where k, is any constant being suitably determined. Next we 
obtain an integral y, of the equation 

al 29 AY = 

dé? ‘sin2@ dé sin?@ ”’ 


(30) 


such that y=0 and dy/d#@=k, at 6=1°, where k is any constant. The conditions 
dy,/d0@~0=~dy,/d6 are very neary satisfied at @= =89° for any solutions. The intervals 
of integration are 1° of colatitude (for @=1°~5°), 2° (for @= 5°~15°), 4° (for @=15° 
~73°), 2° (for 6=73°~81°) and 1° (for 6=81°~89°). Then y=y,+ey is a required 


Fa AND Gn (e.m.u.) 


so 70 860 80 


0 300 40—s«G0 
CO-LATITUDE 
Fig. 3 The variation of f, and gn with cola- Fig. 4 The amplitudes and phases of the elec- 


titude. trostatic field (Xan) 


general integral, and the constant ¢ is determined by the condition (29) at 6=89°. In 
this way the general integrals f, and g, of (10) are obtained as shown in Fig. 3 and 
using (9) and (4) we get the electrostatic field (X, Y) as shown in Fig. 4. In this 
figure and the following figures suffixes 1 and 2 denote diurnal and semi-diurnal com- 
ponents respectively. Substituting these values into (3), we get wind velocities (u, v) 


as shown in Figs. 5 and 6. 
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m/e 


AMPLITUDE 


t) {0 20° 30 50> —§0—"S0> "70> “s0"~" 00 
CO-LATITUDE 


Fig. 5 The amplitudes and phases of the scaloi- 

dal wind system (#n, Un) in the E region, The 

phases denote the local time when the maximum 
values are attained. 
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‘ o. 10: 242— 1 ae een eto 
(2) =r 
( 1000m/SEC (b) 


Fig. 8 Toroidal wind velocities in the E region. 
(a) Diurnal wind system 
(b) Semi-diurnal wind system. 


iz 4 
LOCAL TIME 


t) 10 20 30 40 $0 60 70 80 80 


Fig. 7 The amplitudes and phases of the toroi- 
dal wind system (u,, v,) in the E region. 


6. Application of I.D.A. Analogue Computer 


At first the Integro-Differential-Analyser was used to integrate equations (10) 
and (30) and it was found that a number of solutions exists which satisfy the boundary 
conditions f,=0 (at 6=0) and df,/d@=0 (at 0=7/2). To proceed further as in the 
preceding section, the accuracy of the machine seemed to be a little insufficient, hence 
the integration was carried out as in section 5.* 


7. Drift Velocities in the F Region 


For the present model of atmosphere, the diffusion under gravity is inappreciable 
below 400km level and the vertical drift velocity of ionization is represented approxi- 


mately [12] by 


W=E,cos¢/H (31) 
and the horizontal drift velocities referred to the ground by 
U=—(E,/H)sin¢g+u } (32) 
=(E,/H)(1/sin¢)+v 


The theory of induction drag [3] [5c] [7] indicates that the air in the F region is 
accelerated horizontally by the electrodynamic force due to the electric and magnetic 
fields, in such a way that the dynamo field (vH,, —#H,) cancels most part of the field 
(X, Y). The air is accelerated in a time of order t=p/(o,H,’)=40 min. at the center 


d Mr. A. Kishima of Electrical 
* We are indebted to Professor K. Maeda, Dr. B. Kondo an 
eugadine Institute of Kyoto University, for the use of I.D. A. Analogue computer. When 
the integration was completed, it was noticed that the I.D.A. will be sufficient if the procedure 
of integration is a little altered. But we could not find any chance to repeat the operation. 


% Latitwe 


~ 
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of the F2 region (electron density =10°/c.c.) in the middle latitudes, where o denotes 
the air density. However, the air is little accelerated vertically, since the gravity is 
far effective than the electrodynamic force in this direction. In the present calculation, 


it is presumed that in the F region 


E,=X+vH,=0 } (33) 
E,= Y—uH,=0 
Hence, from (31), we have 
W=—0 (34) 


except near the equator. This means that the value of W is much smaller than that 
produced by electric field Y only. If the semi-diurnal component of E, deviates much 
more from zero than does the diurnal one, then the predominant semi-diurnal variation 
of the #,,F2 in the moderate latitudes at equinox may be interpreted by the drift W 
given by (31). On the magnetic equator, W is directly related to the electric field Y 
and the amplitudes and phases are: 

diurnal 22.8m/sec, ilhs; 

semi-diurnal 4.9m/sec, 12.5 hr; 
and the results are in essential agreement with Hirono and Maeda’s estimation [25]. 
Substituting (33) into (32), we have 


U~u=—(Y/H)\1/sin¢), } (35) 


V~=v=(X/H)(1/sin ¢). 


These values are not much different from 
those which are obtained for #=0=2, 
namely for the F region without air mo- 
tion. The amplitudes and phases of the 
horizontal drift (U, V) are shown in 
figure 9. 


a/m 


8. The Vertical Motion of Ionization in 
the Lower Ionosphere 


In the lower ionosphere, the vertical 4° 
drift velocity of ionization is given [5c] by - 
W=—-7,E,+72ky (36) 
where 20 

71=cot p+(e/m,) {v;/(v2+0,*)}, 
T2=itGe(Di/de), 
p,=(1/m,v,) {v,?cos?$/(v,2 + o,?)+sin?¢}, 


q,=(1/H){o,?/(v,?+,")}cosd, (r= i, e) ow ae go 40 so eo 70 ae 


and —e, v, and w, denote the charge (e.m.u.), cls Coe a phases of the hori- 
collision frequency and the gyrofrequency ocities (Up, Va) in the F2 region. 


of the electron and the suffix i stands for the quantity for the ion. In the F region 


a a en 


a Ps 
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(36) reduces to (31). The equation’ of motion for an electron is dz/dt=W, which may 
be transformed to 


Lay, 
di = 1.37 x 10‘W (37) 


The solution of this equation for a préscribed value of W may be written 
a= '2(2,;'A) (38) 


with 2 fixed by initial conditions at A=4). In this section we confine our attention to 
the region below 250km level during the night. Since in the region, the gas density 
is sufficiently high and the electron density is low during the darkness, the induction 
drag mentioned in the preceding section will be almost negligible. Using the result of 
the case (a) of section 5, we take E,=X and Ey=Y above 160km and in the lower 
region E, and E, given by (28). The curves of (38) are written in Fig. 10 for latitudes 
20° and 45° and for some initial values: 
z=250km and 200km. In these calcula- KM 
tions we did not use the zonal term of 240 
E, which lowers the electrons in higher 220 
latitudes than 35° and lift in the lower 
latitudes except on the equator. More- 
over in the real atmosphere there will be i 180 
the general circulation which appears in o 160 
the expressions of E, and E, as dynamo WwW 
field. It is seen that if the circulation in 140 
the E region is appropriate, some elect- 120 
rons in the F region may descend at 100 
night and enter into the E region. The 

descend on the magnetic equator is al- 

ready discussed in [25]. Recent measure- 

ment [26] by rocket at White Sand 18202224 2 4 6 
(32°20’N, 106°30’W, Geomag. Lat. 41°N) LOCAL TIME 
indicated an existence of N,* ions at Fig. 10 The solutions of the equation dz/dt=W 
midnight between 98km and 120km. Day 0F some ee the night, 
time observation [27] did not show N,* 

in the level. It is suggested that the observed N,* at night may be those which des- 
cended from the F region, where the ionization of N, is likely to take place [18]. 


9. Discussion of the Results 


(1) Electrical Conductivity 
The daily variation (1:12) of the conductivity assumed in section 3 is based on 


the radio observation of the lower ionosphere [28] and on Maeda’s analysis [29], and 
differs from Chapman’s estimation (1:25) based on JZ variation. According to the 
mean S, variation determined by Vestine et al. [30], the ratio of the ranges 4H, 
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(Ohr.~4hr.) to 4H, (12hr.~16hr.) is about 1:13 at many stations in middle latitudes, 
therefore we cannot expect that the night time conductivity is exceedingly small. 
The Z current may originate in a different layer from the S layer. 

Recently Piggot [31] measured the electron density of the E region at night and 
pointed out that the density at midnight is about 1/50 of the noon value. If this ratio 
represented correctly the daily variation of the conductivity of the S bearing layer, 
our present result would have to be greatly revised. The fairly large zonal term Eo, 
obtained in section 4 may be due to incomplete expressions of the conductivity varia- 
tion and the S current system. 

(II) Wind velocity 

Let A, and T,, denote the amplitude of wind velocity in the middle latitudes 
and the time when the wind is directed northward in the northern hemisphere (or 
southward in the southern hemisphere) respectively, and m=1 and 2 denote the diurnal 
and semi-diurnal components respectively. 

(a) Observed Results 

Briggs and Spencer [32] summarise the observed wind by various methods, as 
follows: The mean velocity is about 80m/sec and the periodic part is mainly semi- 
diurnal. In the northern hemisphere A,~30m/sec, 7,~15hr. The rotation of the 
wind velocity is clockwise and anti-clockwise in the northern and southern hemispheres 
respectively. Weiss [33] observed wind in the southern hemisphere using meteor trails 
and concluded that the main periodic wind is diurnal; A,~21m/sec, T,~5hr. 30m. 
and the rotation is anti-clockwise. 

(b) Calculated Results 

For scaloidal wind system 

A,=25m/sec, A,~14m/sec, T,~=7T,~13 hr. 
And for toroidal wind system 

A,=~100m/sec, A,~30m/sec, 7,~7.~13 hr. 
The diurnal wind is much greater than the semi-diurnal one for both the S and T 
winds. The calculated semi-diurnal S velocity is about one half of that observed, and 
the difference of phases between (a) and (b) is to be noticed. The rotation of the S 
wind has the same sense as that observed. For the T wind, the rotation is the same 
as the above wind in the middie latitudes, but opposite in the lower latitudes, where 
the wind in the E region is not yet observed. The observation in the lower latitudes 
will be of importance to decide whether the wind system is nearer to the S or to T 
modes. In the present calculation the condition (33) is presumed in the F region. 
Therefore the daytime effective conductivity of the ionosphere is of order Keatate, 
=7.69 x 10-® except near the equator, to the approximation of [5b]. If there is much 
deviation from (33), for example suppose that ~=0=v in the region above 160 km, then 
the electrostatic field, mainly due to the accumulation of charge by horizontal comno- 
nent of the Hall current, is considerably dissipated through the F region and the 
effective conductivity of the E region takes the value of order a,+4a,2/(a,+ 8,;)=2.5 x 10-8 


in the daytime except near the equator [5b]. Therefore the necessary wind velocity — 


ho 
a 


he SNE tp ahs 
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in the E region for the S, current may have to be enhanced to several times. 

Thus it seems that the calculated wind agrees with the observed values as 
regards order of magnitude. 

(UI) Drift velocity 


In section 7 it is shown that except near the equator, the vertical drift in the 
F2 region is not directly related with the electrostatic field (X, Y) but the horizontal 
drift is determined by the value of (X, Y). Briggs and Spencer [32] show that the 
E-W component in the F2 region is of order 100m/sec and towards the east by day 
and towards the west by night. The calculated horizontal drift velocity (U, V) shown 
in Fig. 9 is mainly diurnal and the phase of V is in phase at high latitudes bu: 
opposite phase at low latitudes with those observed. The velocity is greater than the 
observed velocities in the E region but smaller by a factor 3 than those observed in 
the F2 region. It is to be noticed that in the present approximation, the horizontal 


air motion in the F2 region is approximately equal to the horizontal drift velocity of 
ionization. 


. 


Note: During the course of this investigation, a similar problem was attacked by H. 
Maeda [34] and S. Kato [35] based on an atmospheric model of low density 
determined mainly by Rocket Panel (in which the gas density is of order 108/c.c. 
at 300km). We used an model with higher density (10!°/c.c. at 300km) for strong 
reasons mentioned in section 3, and the methods of solving the dynamo equations 
are essentially different. The wind treated by Maeda is scaloidal. Kato used the 
hydrodynamical equations of motion as a subsidiary condition to the dynamo 
equations. His results concerning the diurnal wind is, however, not yet reliable 
until the mathematical difficulty for singular points of the velocity is overcome. 


10. Concluding Remarks 


From the foregoing analysis the following conclusion may be drown: The 
observed wind velocity in the E region is sufficient to produce the Sy field on the 
ground. The horizontal drift velocity of ionization in the F2 region will be roughly 
the same as the horizontal air velocity. In the middle latitudes, some part of electrons 
in the lower F region will descend at night and enter into the E region when the 
general circulation takes suitable direction. 
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Horizontal Wind Systems in the Ionospheric E region Deduced 
from the Dynamo Theory of the Geomagnetic S, variation 


Part II. Rotating Earth 


a 
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Abstract 


Horizontal wind systems in the ionospheric E region are deduced from the 
dynamo theory by using the same total electric field as in Part I. {1] In Part II 
the effect of the Coriolis force is taken into consideration and a further study is 
made about the wind motion in the E region. The wind velocity is obtained by 
solving the equation of motion of the atmosphere for the rotating earth simul- 
taneously with the dynamo equation. 

It is shown as in Part I that the diurnal wind motion predominates over 
the semi-diurnal one in the E region and the semi-diurnal pressure variation is 
almost in phase with that observed at the ground. 

The conspicious differences in the results between Part I and Part II are: 
(1) The diurnal pressure variation in the E region reaches its maximum value 
at about 7 p.m. local solar time in Part II instead of at 9 a.m. ‘as in Part I. 
(2) Amplification of the semi-diurnal pressure variation in the E region is esti- 
mated to be 25~30 at low latitudes in Part II instead of 50 as in Part I. 

The vertical drift velocity in the F2 region is calculated by using our wind- 


systems. 
In Appendices the mathematical consideration is given to our treatment. 


1. Introduction 


The observation of the daily pressure variation on the ground shows the existence 
of remarkable regularity in the phase and the amplitude of the semi-diurnal variation. 
The theory of the atomospheric oscillation has been studied by G.I. Taylor, C.L. Pekeris, 
K. Weeks, and M.V. Wilkes. According to this theory the atomosphere of the earth has 
a free period of 12 hours [2] [3] [4] [5]. It has been believed that in the ionosphere also 
the semi-diurnal wind motion would be predominant. 

Observation of the wind in the ionosphere has recentely beén conducted by 
using various techniques [6]. However, it is still questionable whether the wind thus 
observed may indicate the real air motion in the ionosphere. 

In Part I H. Maeda deduced the wind system in the E region from the dynamo 
theory of the geomagnetic S, variation. [1] The heterogeneity and the daily-variabili- 
ty of the ionospheric conductivity were taken into consideration in his treatment. He 


( 24) 
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calculated the wind velocity by assuming, for the purpose of simplification, that the 
wind motion is irrotational and associated with a velocity potential. Due to the action 
of the Coriolis force, however, the wind motion is not irrotational. The Coriolis force 
has considerable effect on the diurnal and semi-diurnal wind motions, which play a 
significant role for the geomagnetic S, variation. In Part II, the effect of the Coriolis 
force is taken into account and an advance is made in the treatment of the pro- 
blem. Since no velocity potential exists in this case, the present author makes an 
attempt to obtain the wind velocity by solving the equation of motion of the atomos- 
phere for the rotating earth simultaneously with the dynamo equation. 

Our calculation as in Part I is based on data of the total S, field, internal plus 


external, in the mean solstice, say, : (S+W). 
2. Mathematical Treatment 


The total electric field, corresponding to the geomagnetic S, variation, is denoted 
by E. It consists of the two parts, namely, the electrostatic field and the induced 
electric field, which are indicated as HE, and E; respectively. Then, 


E=E,+E,. (1) 
Operation of curl on both sides of (1) gives 
curl E=curl E,=curl(vx Hy), (e.m.u.) (2) 


because curl E,=0, where v is the wind velocity and H is the geomagnetic field. 
This equation (2) is expressed in the scalar form as follows. 


oe sin 0E,)=  -(@H,) + sin 0uH,), (3) 
where # and v are the southward and eastward components of the wind velocity 
respectively, 8 is co-latitude, and A is longitude. Ey, and Ey, are the southward and 
eastward components of EF, which were calculated by H. Maeda from data of the 
geomagnetic S, variation during the Second Polar Year and the ionospheric conduc- 
tivity by Maeda’s formula. (Part I. Fig. 1A) [1] H, is the vertically upward compo- 
nent of the geomagnetic field, which is assumed to be that of the centered dipole 
alined along the geographical axis. It is 


j= —Ccos8, (4) 


where C is a constant with the value of about * gauss. 
Together with the dynamo equation (3), the equations of motion of the atomos- 
phere for the rotating earth are set as follows. 


— — — ck —— —_—- Ps 

at 2vacosé= — a € ), (5) 
Ov = a 1 0 bp : 

at +2uwcosé=—-—. ‘ant rye ee | (6) 


where is the angular velocity of the earth, @ the radius of the earth, p, the mean 
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static density of the E region, and p’ the perturbing fraction of pressure. The above 
two equations (5) and (6) are readily obtained by assuming that the wind motion is 
produced by a small perturbation” of pressure and the products of small quantities of 
perturbation are neglected in the hydrodynamic equations of motion. 

From these three equations 48), (5), and (6), it is possible to determine D'/ 0 
which corresponds to the geomagnetic S, variation in the ionosphere. If b'/0) which 
is proportional to the pressure variation is known, the wind velocity is to be deduced 


by (5) and (6). 4 
Since the periodic solution is considered here, we take 
BP, erway (7) 


where j= ,/—1, o is an angular frequency, and m is a constant, which takes only 
integers due to the condition of stability. 
Then, from (5) and (6), 


} jo d \E. 
= rr 6 8 
. 4aw* eee do sa Foot Po vee 
me! SE eal sae FE AY 9 
~ 4aw® f2—cos?0 (3 coed dot Sar : (9) 


where f=e/2w. When substitution of (8) and (9) into (3) is made, the right side of 
(3), after some arrangements, becomes 
arg d ial dim 3 oe 
=p eaiCsind| cot Fe Fos sin 0 Ee cae : 
mcos@ 


d m = d,.m 
ee re peal at = ad (10) 
anstebes ~Faeas( 7% = a6 ot zinta) ~ f?—cos*6 do f c 0" 
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The relation between D and 4 is expressed as 
(O-Fn(ovFn)—m(—fjeaw(aen). 00 


- * < © ‘ 
(Dn }Po= Laeger) {vior+(D—Fr ) rot (15) 
From (12) and (14), the righ-hand member of (3), i.e. (10) is equal to 


Oe ceed cence ca a 


The fuctions ¥™, ¥,(™, and ¥,© are assumed to be expanded into a series of the 
associated Legendre functions as 


Then, 


To) = SY Ay m Pas (ity Fa? = 35 Oy m Pay llth, Fa = 31 Boy m Par m(tt), (17) 


where P,, ,(/) is an associated Legendre function of degree m and order m and An, » 
@n,m and By,» are coefficients of expansion. From (11) and (12), the following relations 
are obtained. 


2n-1 Bes — =. 3 


28 : : S. Kato 


; ee Do = : C Ma~ N= 
Left side of (3)= al 00 (sin 0EyY) = sale 2 Fe (Pons Orne SETA). (21) : 
Right side of (3)= 0 (vH.) + 0 (sin ; nH any ‘\ vy P pjm(we+ dr) 
Ri Maen OF ean Veen ee a of! 
Fsas % 2(7 | 
where Fob as m*sin@ .| n—m 29m _y ee leer 3 
2 afw 2n—] f m Pn-1y)m™ 
n(n+1)— f 
‘mr * 
9 
n+m+1]., ae 2( *) . 
: n2+2n4+2™ 4 \S po (23) 
2n+3 F m fact te 
ate 
Table 1 The spherical harmonic coefficients of (21) 
Fy,m is expressed as Fy.q=F\y+jF Oy 
(unit: 10!°cm*sec™') 
M N | FW | Fou 
1 2 189.1 —141.3 
4 48.08 —13.33 | 
6 10.26 — 40.41 
8 13.65 —12.97 
2 3 — 22.74 11.98 
5 —3.733 —2.414 
7 3.024 —2.701 
9 2.906 — 1.982 


and the value of Fy,» is given in Table J. In the above equation (22), c=mw is 

assumed, because the S, variation is dependent on local solar time. The appearance 

of 7 before the bracket of (23) shows phase difterence of 90° between the both sides of (3). 
If w tends to zero, it is readily shown from (20) and (23) that 


a 


me ieee | (n—m)(n+1)(n—)) 4 (ntm+l)n(n+2) 4 | 

.. inte 5 — eee Cie m—ly m : = n : me . 24 
+ = 2n—1 im f 2n+3 tirm |. (24) 
This is just identical with the expression in case of non-rotating eath, treated in Part 


I, in which it is assumed that 
v=—grad 0( = : grad a (25) 


00 

where @ is the velocity potential. [1] From (21) and (22), Bn, » is determined by equat- 
ing the coefficients of each P,, ,,() and substitution of the values of 7 into (20) 
makes it possible to evaluate A,,,,, the amplitude and the phase of pressure variations 
required. In the actual calculation, the first 4 terms of Fy, mw and L,,,, are considered. 

The soulution obtained here is a particular integral of the differential equation 
(3). Therefore, it is not yet certain that this solution is unique. The uniquness of the 
solution we have obtained will be proved in Appendix 1. In Appendix II the mathema- 
tical consideration is given to the singular points of the dynamo differential equation, 
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3. Pressure Variations and Wind Systems 
The amplitudes and phases of various harmonic components of the pressure 
variations are shown in Table 2. The pressure variations are expressed in the form as 


Table 2. The spherical harmonic coefficients Pes . =STAn, mPn, mSin (Mut + Mm). + an, m) 


(Unit of An,m: 10°cm*sec™? Unit of an,m: degree) 
m n Aase a 
1 | 1 114.1 150 30’ 
| 3 33.2 | 141° 
| > 3.30 200° 
2 | 2 8.79 135°30 
4 0.567 | 288° 
| 6 0.00334 152°45 
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ee seo time na Local solar time (hours) 
Fig. 1 Distribution of p’/p, which is propor- Fig. 2 Distribution of the wind velocity corres- 
tional to the pressure variation. Diurnal com- ponding to the pressure variation as shown in 
ponent. Unit: 10°cm?sec~’. Fig. 1. Arrow-scale: 50m/sec for a length 


equal to the side of the small squares. 
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Local solar time (hours) Local solar time (hours) 
Fig. 3. Distribution of p'/p, which is propor- Fig. 4 Distribution of the wind velocity corres- 
tional to the pressure variation. Semi-diurnal ponding to the pressure variation as shown in 
component. Unit: 10°cm?’sec™?. Fig. 3. Arrow-scale: 50m/sec for a length — 
equal to the side of the small squares. 
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fe, ph time (hours) Local solar time (hours) 
Fig. 5 Distribution of p’/p, which is propor- Fig. 6 Distribution of the wind velocity corres- 
tional to the pressure variation. Sum of the ponding to the pressure variation as shown in 
diurnal and semi-diurnal components. Unit: Fig. 5. Arrow-scale: 50m/sec for a length 
10®cm?sec™?. equal to the side of the smal squares. 


se 
ssa 
~~ 
PASS 
= 


2_yA,,. Pa, m(z)Sin (mwt+mA+an, m), in Which A,,,, is the amplitudes and a,,,, the 
phase angles. The global distribution of the pressure variations (isobars) and the wind 
velocity corresponding to those pressure variations are given in Figs. 1 to 6. The 
figures show the distributions in the northern hemisphere. The distributions in the 
northern and southern hemispheres are symmetrical about the equator. 

The remarkable results of our calculation are: 

(1) The predominant component in the pressure variations in the E region is 
not semi-diurnal, but diurnal. The diurnal component is about 3.5 times as great as 
the semi-diurnal one. 

(2) The semi-diurnal pressure variation in the E region is almost in phase with 
that observed at the ground. The maximum occurs at 10 am. and 10 p.m. local 


solar time. 
(3) The diurnal pressure variation reaches the maximum value at 7 p.m. local 


solar time. 
(4) The amplification of the pressure variations is estimated by using the follow- 


(5) 

Do/i 
(fe 
Po/s 

where pf, is the static pressure and T the static temperature. The subscripts 7 and g 

refer to the ionosphere and the ground respectively. The mean molecular weight in 


- the ionosphere and that near the ground are assumed to be equal. Provided 7,—T,, 
p’~0.937 mmHg, [5] and p.=0.001293 g cm~* near the ground, we get 


ing relation. 


(Fh 
Af0"*t ~25~ 30, 
(5), 
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for the semi-diurnal case at low latitudes. 

(5) The wind velocity is 7~8m/sec. for the semi-diurnal case and 15 
~20 m/sec. for the diurnal case at low latitudes. At high latitudes the wind velocity 
becomes greater for both cases. 

The first two of the results mentioned above ‘are almost the same as in Part I, 
in which the rotation of the earth is not taken into account. However, the phase of 
the pressure variation and the wind velocity in the diurnal case are far different from 
those obtained in Part II. [1] (See Part I. Fig. 2B) 

Such a situation seems to be acceptable because the effect of the rotation of the 
earth becomes greater with the increase of period of the pressure variation. In the 
case of non-rotating earth, the maximum of the diurnal pressure variation occurs at 
9 a.m. instead of at 7 p.m. asin Part II. The wind velocities in Part I are about twice 
as great as those in Part II for both diurnal and semi-diurnal cases at low latitudes. 
The diurnal wind system in Part II is cyclonic in the low pressure zone and anti- 
cyclonic in the high pressure zone, in the northern hemisphere. 

M. Hirono and T. Kitamura deduced the wind system in the ionospheric E 
region from the dynamo theory. [8] They used the similar data of the geomagne- 
tic S, variation as we, but made some different assumptions about the geomagnetic 
S, variation at the pole and the conductivity at night. They considered the following 
two cases separately. (1) The wind is entirely scaloidal. (2) The wind is entirely 
toroidal. The actual wind, however, must be of intermediate character between these 
two extreme cases. It will be noted that this is the case in Part II, where the 
Coriolis force may be a cause for toroidal motion. For instance, in the diurnal case 
at high latitudes, where the Coriolis force is very effective, the toroidal wind becomes 
predominat, while at low latitudes both toroidal and scaloidal winds are approximately 
equal in magnitude and opposite in sense. Hence, at low latitudes the winds of these 
two kinds nearly cancel out, and the small toroidal part is left. In the semi-diurnal 
case the Coriolis force is not so effective as in the diurnal, and accordingly the scaloidal 
part becomes greater than the toroidal one. However, the difference of the semi-diurnal 
wind velocity in magnitude between Part I and Part II suggests the existence of 
considerable contribution of the toroidal part to the wind motion. 


4, Concluding Remarks 


The results of the recent measurements of the horizontal wind velocity in the 
ionosphere were reported by B.H. Briggs and M. Spencer. [9] The semi-diurnal wind 
obtained from our calculation agrees with those results. In order to test our calcula- 


tion the observation of the wind in the ionosphere-must be conducted at various 


latitudes. 

As an application of our calculation it seems worthwhile to estimate the 
electrodynamic drift velocity in the vertical direction in the F2 region. The vertical 
drift velocity W is given as [10] 


a 
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iWw= = cos¢, 


where F is the intensity of thé geomagnetic field, ¢ the dip angle of the geo- 
magnetic field, and E, the eastward component of the electric field in the F2 region. 
We assume that there is no wind motion in the F2 region. Then, the electric field in 
the F2 region becomes electrostatic. Since the electrostatic field in the F2 region is 
approximately equal to that in the E region, it is estimated from our calculation. 
Substraction of the induced field from the total -electric field gives the electrostatic 
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Fig. 7 Distribution of the diurnal static field. Arrow-scale: 10%e.m.u. 
for_a,length7equal to the side of the small squares. 
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Fig. 8 Distribution of the semi-diurnal static field. Arrow-scale 
10°e.m.u. for a length equal to the side of the small squares. 


-10 


Fig. 9 Daily variations of vertical drift 
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field. The global distribution of the electro- 
static field is shown in Figs. 7 and 8. The 
daily variation of the drift veocity at @=50°, 
70°, 80°, and 90° is given in Fig. 9. 
intensity of the geomagnetic field at the 
From this 


with such a character. 
12 14 #16 18° 20 22° 24 


We have assumed that the ionosphere 
And the treatment 
is two-dimensional. Therefore, the pressure 


is a thin spherical shell. 


variation and the wind system in our case 

should show some average state of the entire 
ionosphere, though much weight is put in the 

20 E region, on account of its greater conduc- 
tivity. 

10 : 
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region makes the phase of the atomospheric 


concerned. 


Local Solar: Time 


action in the ionosphere. 


velocity of electrons in the F2 region at : : 


6=50°, 


70°, 


80° and 90°. 


equator is taken to be 0.28 gauss. 
figure 9 it is known that the vertical drift 
velocity is principally diurnal and the time of 
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occurrence of the maximum velocity shifts 
from 11 a.m. at the equator to 8 a.m. at 0< 
70°. The geomagnetic distortion of the F2 
region near the magnetic equator will well 
be explained by the vertical drift motion 


Actually, the existence of the positive 
temperature gradient at the bottom of the E 


oscillation change rapidly with height. [8] [9] 
This actual situation requires a three-dimen- 
sional treatment. However, the outstanding 
features of our results will remain valid, 
even if more complicated treatment is made, 
as far as the geomagnetic S, variation is 


The diurnal atmospheric oscillation has 
never been studied as yet. Our resuts seem 
2 44 16 YS 20 2 2 to Suggest that the diurnal wind is playing 
an important role in the ionosphere. 
required that a study of the diurnal oscilla- 
tion be pursued in future. The diurnal wind 
motion may be produced by some thermal 
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APPENDIX I 


Uniquness of our solution of the dynamo differential equation 


The dynamo equation has been solved by expanding both sides of the equation 
into a series of the associated Legendre functions and equating the coefficients of each 
harmonics. 

The solution v, thus obtained is a particular integral of the differential equation. 
The general solution is 

v=0,+kv,, (A1.1) 
where v, is the general solution of the homogeneous equation, which is obtained by 
putting the right side of the equation to zero, viz. 


0 —_ ; 
: a = Al. 
aa (v.H,) + a0 (sin 6u,H,)=0, (A1.2) 


and & is an integration constant to be determined by a boundary condition. 

What follows will prove that the homogeneous equation has no solution to 
satisfy our boundary condition, and, therefore, the particular solution, i.e. the solution 
obtained in our treatment is unique. 

The necessary and sufficient condition to be satisfied by v, in (A1.2) is 


1 OA 


—uwH, Se Ss A)? 
pe sds (A1.3) 
0,1, a 90” 


where A is some scalar function, which is analytic in the domain of our discussion 
and is arbitrary in this equation. H, is assumed to be 
* 


H,=—Ccos8, 
as (4) in section 2. Then, 
- 1 0A 
“~~ CcosO@sind OA’ nee 
gece 1A 
>" Ccosé 00° 


These are substituted into (5) and (6) in section 2. After elimination of the terms 
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of pressure, (5) and (6) give the following relation. 
0 f Ou De Ov Oy i | 
—}——— — sin Ov) }—2 7] sin 6 cos 62) }=0 Al.4 
a ee wy) (sin »} ofcos aA + a0 (sin 6 cc wu) ( ) 


Since only periodic solutions are to be considered, substitution of (A1.3) into (A1.4) 
gives 
aA dA m? (AL.5) 


tan 6 — A=0 
doz . do sin?0 


Here A « e4(%+™) is assumed as (7) in section 2. Next, transformation is made like 


z=sin??, A=22w. 


Then, (5) becomes 


(2-1) 4 + m™+1)2-c1+m) |e + (F)(5 )e=0- (A1.6) 


The above equation is known as the hypergeometric equation, the solution of which 


is expressed as follows. [11] 
w=(%, 7 1+ m| 2), (A1.7) 


This function converges in the closed domain 0<z<1. At z=1, that is, at the equa- 


tor, [11] ra ; 
x +m 


where 7’ stands for the gamma function. From the first equation of (A1.3’) 


(A1.8) 


im 
cu,cosé= — =A (A1.9) 


The relation of (A1.9) requires that A=0 at 0=90°, provided z, is finite at the equa-. 


tor. This boundary condition is satisfied only by R=0 in (Al1.1). 

It is concluded that the solution obtained in section 2 is an unique solution, 
which is physically available. It is equivalent to saying that the induced electric field 
has no scalar potential. The wind velocity is determined from the dynamo theory 
without any mathematical arbitrariness. 


APPENDIX II 
The wind velocity about singular points of the dynamo differential equation 


The differential equation (3) in section 2 contains three singular points which 
are located at cos@=0, f and 1. As our mathematical treatment shows in section 2, 
we have obtained the analytic solution which is expressed in a series of the associated 
Legendre functions. Therefore, the solution, i.e. the pressure variation, contains no 
singular point. However, if we substitute the pressure variation obtained into (8) and 
(9) in order to get the wind velocity, the wind velocity becomes indeterminate or 


infinitely great about the singular points mentioned above. Such a fatal situation never 
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occurs if we use the auxiliary functions, ¥Y, and ¥, in section 2, for evaluation of the 
wind velocity. 


What follows will show the-method by which it is possible to determine the 
wind velocity about the singular points. 


For the purpose of simplifiation, we put the wind velocity as follows. 


U=U'(A)eiettmr) | (A2.1) 
v=v/(A)eiottmr) (A2.2) 


where zw and wv’ are the function of 8. Then, from (8) and (9) in section 2, (A2.1) 


becomes 
as - { ss uae cot o}rom 


~ haw? f*—cos’6 | dé 


eid u (p-™ ae A2.3) 
4 aw* sin 0(f?—cos?0) fi . 
_ —jo 1 a) (D-2 \y Pe 
haw siné f? | . rig : 


We expand ¥, and ¥, into a series of the associated Legendre functions as (17) in 
section 2. Then, 


ie : 
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“(n+ 1) 
_ (nt mt V(ntm+2) oF pa ee ose, P,, (tt) ‘i 
(2n+3)(2n+5) = Ee m(n+3) lava =e 6° (A2.6) 


In (A2.4) and (A2.6) the denominators never vanish in our calculation, in which 
Bits Baits Bist s 2+ *49 Buty Pers Bea .o>ss, are COnsidered.< farthec, 
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6+0 sin@ 
\ 


= finite. 


From the above argument, we can know that it is possible to evaluate the wind velo- 
city at any point between 6=0° and 90°. 
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